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Abstract. In the best observed spiral galaxies, the HI is observed to 
decline slowly with radius for up to 90% of the maximum extent, before 
experiencing a rapid truncation. Early models predicted that the outer 
parts of spirals would be fully ionized by a metagalactic UV field at 
column depths of roughly 10^^ atoms cm~^. Weak optical line emission 
has now been observed beyond the HI edge in a number of cases and, 
furthermore, this warm gas has been used to trace the dark halo far 
beyond the HI cut-off. But the observed emission measures (40—90 mR) 
are higher than expected for a cosmic origin. A more likely explanation 
is ionization of the warped outer disk by the blue stellar population at 
smaller radius. There are many implications of ionized warped edges, in 
particular, severe HI warps may be greatly under represented in current 
HI surveys. We briefly discuss the implications of this model for the 
Galactic warp. 



1. Introduction 

Bochkarev & Sunyaev (1977) first argued that the HI disks of spiral galaxies 
truncate at a few times the optical diameter because the exponentially declining 
HI column eventually becomes fully ionized by a metagalactic UV field. After 
detailed calculations of the expected emission measures (Maloney 1993; Dove 
& Shull 1994), it was soon realized that these extremely faint levels could be 
reached with the Fabry-Perot 'staring' technique (Bland-Hawthorn et al. 1994; 
hereafter BTVS). This method has since been used successfully on Sculptor 
galaxies (Bland-Hawthorn, Freeman & Quinn 1997; hereafter BFQ), M33 and 
NGC 3198 (Bland-Hawthorn, Veilleux & Carignan 1998; hereafter BVC). In 
keeping with the main focus of this workshop, we start with a brief discussion of 
the dark halo in NGC 253. But most of the article will be devoted to the source 
of the HI truncation. 



2. The dark halo of NGC 253 

In Fig. la—c, we show the kinematic measurements along the major axis of NGC 
253 deduced from VLA HI (inner points), TAURUS-2 Ha and [Nil] lines (outer 
4 points). The HI data are from a fit to the full velocity field (Puche, Carignan 
& van Gorkom 1991; hereafter PCvG). Whether one fits to the approaching 
side, receding side, or the full velocity field, the rotation curve is still rising in 
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Figure 1. In the left panel (a — c), the HI, Ha and [Nil] rotation 
curves are shown with representative models for NGC 253 (see text). 
In the right panel (d — f), we illustrate the influence of an ionizing disk 
on different warp extremes in M33. We show the projected HI, seen 
from the galaxy nucleus, before (d) and after (e) the disk ionization 
is switched on, and the resulting trend (/) in emission measure with 
radius. The deprojected HI data are from Corbelli &: Schneider (1997). 

the outermost parts, a phenomenon observed in a number of low mass spirals 
(Puchc & Carignan 1991). 

The inflexion in the rotation curve is understood naturally to arise from 
the potential of an exponential stellar disk with a much larger dark halo (e.g., 
Carignan &: Freeman 1985). To illustrate this (Fig. la — c), we show represen- 
tative fits to the PCvG and TAURUS-2 data using a 3-component mass model 
(disk, halo, gas). The contribution from the HI surface density (dotted line) is 
the same in all models. In a, we adopt an exponential disk (short-dashed line) 
and choose an isothermal sphere for the spherical halo (long-dashed line) using 
the numerical approach developed by Carignan (1985). In b, we have used the 
disk-halo models of Dehnen & Binney (1997). In c, we have truncated the halo 
at a radius of 10 kpc to demonstrate a possible explanation for the last measured 
point. (A detailed discussion is given in BFQ.) 

The outermost TAURUS-2 Ha measurement suggests that the rotation 
curve may be falling beyond a radius of about 10 kpc, but there are certainly 
other possible explanations for its low observed velocity (e.g., tidal forces). But 
BFQ consider the possibility that the rotation curve of NGC 253 is indeed falling 
beyond 10 kpc due to a truncated dark halo (see Fig. Ic). This compares to 
the much larger dark halo distributions observed in our Galaxy (e.g., Freeman 
1996) and other large spirals (e.g., Zaritsky et al. 1996). 

If this interpretation is correct, there are some interesting consequences. 
The properties of dark halos are best studied in disk galaxies for which the HI 
distribution extends well beyond the optical distribution. In the outer regions 
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of these galaxies, the ratio of (dark matter surface density) / (HI surface density) 
is roughly constant {e.g., Bosma 1978; Carignan 1991). NGC 253 is not such 
a galaxy. Its HI extends only 40% further than the optical extent. But the 
data provide the first hint that where the HI and the light are co-extensive, the 
dark matter also may not extend much beyond the optical distribution. It is 
tentative evidence of the apparent link between the dark matter and the HI. 
We emphasize that the dark matter is still essential to generate the observed 
rotation curve for NGC 253: the inferred mass ratio of dark matter to luminous 
matter is about 5. 

3. What ionizes the outer edge of NGC 253? 

In NGC 253, we see ionized gas with ~ 40—90 mR at and beyond the HI 
edge. The [NII]A6548/Ha ratio is close to unity which is anomalously high for 
diffuse gas. 

Young stellar disk. BFQ attempt to explain this enhancement with recourse 
to standard ionization codes and an idealised model for the young disk (see 
Bland-Hawthorn &: Maloney 1997; 1998; hereafter BM). A useful approxima- 
tion to this model, and the predicted halo fluxes, are given in the appendix. 
Since essentially all spirals are warped, it is likely that this interpretation has 
application to a large fraction of objects. P.R. Maloney has calculated the HI 
truncation produced by our model in the context of M33 (Fig. Id — /). The 
curves correspond to integral sign warps with a range of maximum amplitudes: 
10° (solid), 20° (dot), 40° (short dash), 60° (long dash), modified 60° (dot-dash). 
The last curve is the same as the previous curve except that the column density 
has been reduced by 5%. The existence of a truncated edge is very sensitive to 
the ionizing flux and the shape of the warp. 

An important prediction of the models is that most of the low ionization 
lines are greatly enhanced with respect to Ha, and that [OIII]5007 is essentially 
invisible. This has largely been confirmed from a run in October 1997 with 
K.C. Freeman at the MSO 2.3m using the Double Beam Spectrograph. We 
obtained a total exposure time of 5.6 hours at field positions SW3 and SW4 
(BFQ). The stacked data reveal that, like [NII]/Ha, [SII]/Ha is extraordinarily 
high (~ 0.5 — 0.8) in agreement with the model. 

Cosmic ionizing background. From Ha non-detections towards extragalactic 
HI, the current best 2a upper limit (Vogel et al. 1995) for the cosmic ionizing 
background is J°2i < 0.08 ($° <2-4xl0^ phot cm~^ s~^; see appendix). This 
compares with J°2i ~ 0.006 deduced from the proximity effect (Kulkarni & Fall 
1995). If the 3C 273 sight line is truly representative, the Milky Way dominates 
the ionizing field out to 400 kpc on average, and maybe 700 kpc along the polar 
axis. This is an order of magnitude smaller than the scale of typical galaxy 
separations at the present epoch. 

BTVS show that the expected emission measure from a cloud is unlikely to 
be detectable beyond 300 kpc, or maybe 500 kpc along the polar axis. This has a 
number of important ramifications, e.g., it is unlikely we can detect the mutual 
ionization of M31 and the Galaxy. At this meeting, L. Blitz has laid claim to 
a large fraction of high velocity clouds, placing them at extragalactic distances 
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Figure 2. Ionized outer regions arising from (a) a simple warp, (b) a 
precessing warp, and (c) a simple warp that folds back on itself. 



within the Local Group. However, if most of these clouds produce detectable 
Ha, this interpretation is unlikely to be correct. 

Suprathermal particle heating. The enhanced low ionization lines require 
substantial heating of the gas without further ionization. BFQ go some way 
to achieving this with a dilute, hardened radiation field impinging on gas with 
refractory element depletion. But the required heating rate is too low by at least 
a factor of two. Many early candidates {e.g., relativistic particles) have proved 
unsuccessful. More promising are subrelativistic, heavy (suprathermal) particles 
channelled along magnetic field lines to the outer disk. Their stopping lengths 
can be very long so that the ionization rate is kept low, and a single nucleus can 
produce many energetic "knock on" electrons. 



4. HI disk morphology 

There are a number of interesting consequences arising from global ionization 
of the outer disk. Here, we illustrate how different sources could produce quite 
distinct behaviour in (i) the outer HI disk morphology, and (ii) the trend in 
with radius. In the next section, we consider the special case of the Galactic 
warp. 

In the table below, the dcprojectcd shape of the HI disk at a fixed column 
density is shown in the first column. The source of the truncation and the 
expected trend in Sm (with increasing radius) are given in the second and third 
columns respectively: 



shape of HI disk source of truncation radial trend 

mirror symmetric integral-sign warp + internal source £m increasing 

axisymmetric ambient radiation field decreasing 

bisymmetric tides due to external galaxy £^ = Q 

asymmetric ram pressure from external medium £m complex 
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The first interpretation assumes a simple warp for the HI, a characteristic 
of most warps (Briggs 1990), although one can certainly construct different 
scenarios (see Fig. 2) for ionization by a central source. The prediction of 
£„i increasing is confirmed in Fig. 1/ although only in the limit of strong warps 
(> 30°). These models predict azimuthal incompleteness in the HI distribution 
at critical velocities. Careful kinematic analysis of the best observed spirals 
could put strong limits on each of these models. The second interpretation is 
the conventional picture in terms of a mctagalactic UV field (sec §1) but, to 
date, we see no evidence for such a model. The third and fourth models are 
speculative at best: it is difficult to see how either could produce well defined 
HI truncations. 

In order to explore these different cases, we are currently studying a range 
of galaxy types in different environments, with the aim of tracing the dark 
matter beyond the HI limit. The galaxies chosen for detailed study include 
M31, M33, M81, M83, NGC 628, NGC 3198, NGC 5266, Fourcade-Figueroa, 
and the Sculptor Group. The most detailed work has concentrated on NGC 253 
(BFQ), M33 and NGC 3198 (BVC). 

5. The Galaxy 

In Fig. 3, we illustrate the predicted halo ionizing field. The 21 cm data are 
from Burton (1984) and Burton & te Lintel Hekkert (1986). This is normalized 
to explain the Magellanic Ha stream detections (Weiner & Williams 1996). The 
northern HI warp is easily truncated by the disk radiation; we suggest that the 
southern warp is much less pronounced due to the ionizing influence of both the 
LMC and the inner disk. 

The predicted levels of Ha emission should be easily detectable (150—200 
mR) with WHAM, and only marginally by other all-sky Ha surveys. We antic- 
ipate that the projected Ha emission off the plane will show peaks at I = 70° 
and / = 250° respectively, with the southern extension being more pronounced. 
It should be possible to extend the Galactic rotation curve beyond the current 
HI limit. More generally, it seems highly likely that severe warps in spirals, 
particularly those with UV-bright inner disks, are suppressed by disk ionization. 

6. Discussion 

The ionization models for warped edges suggest that essentially any spiral with 
a young, stellar population can be reasonably traced to and beyond the HI limit. 
In the months to come, there will be half a dozen spirals with optical rotation 
curves at and beyond the HI edge. For these objects, the total mass in H+ is 
comparable to the HI mass. 

It is possible to use the 'staring' method down to surface brightness levels at 
least 6 magnitudes below sky in a lA band (BTVS).It may be possible to extend 
this limit further with telescope nodding, charge shuffling and deep depletion 
CCD arrays, all of which will be in place at the AAT by early 1998. In principle, 
stellar rotation curves can be traced to several magnitudes below sky, although 
the Lorentzian instrumental profile is a handicap to reaching much beyond a 
Holmberg radius! There are good reasons for pushing the stellar rotation curves. 
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Figure 3. The galactic halo ionizing field. The coordinates are with 
respect to a plane perpendicular to the Galactic disk (with the Galactic 
Centre at the origin) at a constant galactic azimuth angle (75°, 255°). 
The dotted lines show the ionizing flux (1^4) due to the galactic disk; 
the solid lines include the contribution from the LMC. The opacity 
of the HI disk (shown in half tone) has been included. The contours, 
from outside in, are for logip4 = 1,1.25,1.5,1.75,2,2.25,2.5,3 phot 
cm~^ s~^. The minor contribution from the Galactic corona is omitted 
(BM). 
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In special cases, the H"*" kinematics may need second order corrections {e.g., 
asymmetric drift) for ionization effects {e.g., Yorke, Bodenheimer &; Tenorio- 
Tagle 1982; Yorke & Welz 1996). 

We stress that 21 cm mapping of spirals is crucial to this work, and is sub- 
stantially more efficient than 'staring' over the HI extent. (For reasons outlined 
in §4, it is even worth the effort to revisit archived 21 cm data.) It takes 0.5—1 
nights to obtain a pair or triplet of optical kinematic measurements. More rapid 
progress is possible with a dedicated 2.5m— 5m telescope (see BFQ, eqn. 7). 
These measurements are tied to the HI data and must be corrected for warping 
(and flaring) of the disk. While Fabry-Perot 'staring' has extraordinary potential 
for many problems, e.g., R.J. Reynolds' outstanding work with the WHAM in- 
strument, it merely complements the deepest measurements obtainable in other 
bands. 

Acknowledgments. 1 am indebted to K.C. Freeman, P.R. Maloney, H.W. 
Yorke, W.B. Burton, and M.A. Dopita for assistance with several aspects of this 
work. 

References 

Bland-Hawthorn, J., Taylor, K., Veilleux, S. & ShopbeU, P.L. 1994, ApJ, 437, 
L95 

Bland-Hawthorn, J., Freeman, K.C. & Quinn, P.J. 1997, ApJ, 490, 143 

Bland-Hawthorn, J. & Maloney, P.R. 1997, PASA, 14, 59 

Bland-Hawthorn, J. & Maloney, P.R. 1998, ApJ, submitted 

Bland-Hawthorn, J., Veilleux, S. &; Carignan, C. 1998, ApJ, in preparation 

Bochkarev, N.G. k Sunyaev, R.A. 1977, Soviet Astr., 21, 542 (originally 1975, 
AZh, 54, 957) 

Bosma, A. 1978, PhD, University of Groningen 

Briggs, F. 1990, ApJ, 352, 15 

Burton, W.B. 1984, in Galactic & Extragalactic Radio Astronomy, eds. K.I. 

Kellerman and G.L. Verschuur (Springer- Ver lag), p. 82 
Burton, W.B. & tc Lintel Hekkert, P. 1986, A&AS, 65, 427 

Carignan, C. 1985, ApJ, 299, 59 

Carignan, C. 1991, contribution to Space Telescope Workshop on Dark Matter. 
Carignan, C. & Freeman, K.C. 1985, ApJ, 294, 494 
Corbelh, E. & Schneider, S.E. 1997, ApJ, 479, 244 
Dehnen, W. & Binney, J. 1997, MNRAS, 287, 5 
Dove, J.B. & Shull, J.M. 1994, ApJ, 423, 196 

Freeman, K.C. 1996, in Unsolved Problems of the Milky Way, ed L. Blitz and 

P. Teuben (Dordrecht: Kluwer), p. 645. 
Kulkarni, V.P. & Fall, M. 1995, ApJ, 453, 65 
Maloney, P.R. 1993, ApJ, 414, 41 
Puche, D. & Carignan, C. 1991, ApJ, 378, 487 
Puche, D., Carignan, C. k van Gorkom, J. 1991, AJ, 101, 456 

7 



Vogel, S.N., Weymann, R., Rauch, M. k Hamilton, T. 1995, ApJ, 441, 162 
Weiner, B.J. & Williams, T.B. 1996, AJ, 111, 1156 
Yorkc, H., Bodcnhcimcr, P. & Tcnorio-Tagle, G. 1982, A&A, 108, 25 
Yorke, H. & Welz, A. 1996, A&A, 315, 555 

Zaritsky, D., Smith, R., Prenk, C, & White, S. 1996, ApJ, 478, 39 



A Ionization models 



We present a useful approximation to the Galactic halo ionizing field deduced 
by BM. We assume an electron temperature Tg ^ lO^K, as expected for gas 
photoionized by stellar sources, for which the Case B hydrogen recombina- 
tion coefficient is ^ 2.6 x 10~^^(10^/re)°-^^ cm^ s~^. At these tempera- 
tures, collisional ionization processes are negligible. In this case, the column 
recombination rate in equilibrium must equal the normally incident ionizing 
photon flux, a^neNfj+ = (fi, where is the rate at which Lyc photons ar- 
rive at the cloud surface (photons cm^^ s~^), is the electron density and 
Nff+ is the column density of ionized hydrogen. The emission measure is 
£-m = J n(,nfj+ dl = Uf^riu+L cm^^ pc where L is the thickness of the ion- 
ized region. The resulting emission measure for an ionizing flux ipi is then 
£m = 1-25 X 10^^(/74 cm^^ pc (= A.5(p4 mR) where (pi = 10^(yO4. For an optically 
thin cloud in an isotropic radiation field, the solid angle from which radiation 
is received is O = 47r, while for one-sided illumination, $7 = 27r. For our disk 
model, however, is anisotropic and can be considerably less than 2it. 

For the cosmic field, J°2i the ionizing flux density of the cosmic back- 
ground at the Lyman limit in units of 10~^^ erg cm~^ s~^ Hz~^ sr~^; 
(= 7rJ°2i/^) is the equivalent photon flux at face of a uniform, optically thick 
slab. 

For the Galactic halo (Fig. 3), to a good approximation, the poloidal UV 
radiation fleld is 

(P4 = 2.8 X 10^ e-^ rkpe-2 cos°-^^+°-^ phot cm-^ s'^ (1) 

where Q is the polar angle (0 < 6 < |) and the Lyman limit optical depth 
T < 10. It follows that the solid-angle averaged flux is 

2.8 X 10^ e"^ , _2 _i 

^"^ = (0.6r + 1.5)r,pe^ ^" " 

In order to explain Magellanic Stream Ha detections, our model favours r « 2.8. 
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